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The g'g relativistic interaction in the Wilson loop approach 
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We study the qq relativistic interaction starting from the Feynman-Schwinger representation of the gauge- 
invariant quark-antiquark Green function. We focus on the one-body limit and discuss the obtained non- 
perturbative interaction kernel of the Dirac equation. 



1. INTRODUCTION 

The dynamics of a system composed by two 
heavy quarks is well understood in terms of a po- 
tential interaction (static plus relativistic correc- 
tions) obtained from the scmirclativistic reduc- 
tion of the QCD dynamics 

If at least one of the quarks is light the sys- 
tem behaves relativistically (the light quark can 
no more be considered static) and a pure rela- 
tivistic treatment becomes necessary (via Dirac or 
Bethe-Salpeter equations). A lot of phenomeno- 
logically justified relativistic equations can be 
found in the literature but up to now we miss a 
relativistic treatment which follows directly from 
QCD. Our work goes in this direction. In order 
to simplify the problem we focus on the heavy- 
light mesons in the non-recoil limit (i. e. in- 
finitely heavy antiquark). Only at the end we 
will briefly discuss the two-body case. Our start- 
ing point is the quark-antiquark gauge-invariant 
Green function taken in the infinite mass limit 
of one particle. The only dynamical assumption 
is on the behaviour of the Wilson loop (i. e. on 
the nature of the non-perturbative vacuum). The 
gauge invariance of the formalism guarantees that 
the relevant physical information are preserved at 
any step of our derivation. In this way we obtain 
a QCD justified fully relativistic interaction ker- 
nel for the quark in the infinite mass limit of the 
antiquark. This kernel reduces in some region of 
the physical parameters to the heavy quark mass 
potential, and leads in some other region to the 
heavy quark sum rules results, providing in this 



way an unified description. We discuss our result 
with respect to the old-standing problem of the 
Lorentz structure of the Dirac kernel for a confin- 
ing interaction 1^-^. The main results presented 
here can be found in 

2. THE RELATIVISTIC INTERACTION 
IN THE ONE-BODY LIMIT 

The quark-antiquark Green function is given in 
the quenched approximation by 



Ginv{x,u,y,v) = (Tr i S^^\x,y;A)U{y,v) 
X iS^'^\v,u;A)U{u,x)) , (1) 



where the points x,y,u,v are defined as in Fig. 
^ ( ) means the normalized average over the 
gauge field A^, S^'^ is the fermion propagator in 
the external field A^, associated with the particle 
i and the strings 

U{y,x) = 

P exp l^ig J ds{y- x)''A^{x + s{y - x))^ , 

are needed in order to have gauge invariant initial 
and final bound states. A very convenient way to 
represent it is the so-called Feynman-Schwinger 
representation (see and refs. therein), 

where the fermion propagators are expressed in 
terms of quantomechanical path integrals over the 
quark trajectories (zi(ii) and ^2(^2)) 
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From Eq. (||) it emerges quite manifestly that 
the entire dynamics of the system depends on the 
Wilson loop: 



(3) 



W{T;A) = TrP expjig (j) , 

being F the closed curve defined by the quark 
trajectories and the endpoint strings U{y,v) and 
?7(m, x). 

In order to treat a simpler case, let us as- 
sume, now, that the antiquark moving on the sec- 
ond fermion line becomes infinitely heavy. The 
only trajectory surviving in the path integral of 
Eq. (||) associated with the second particle is the 
static straight line propagating from v to u. The 
corresponding Wilson loop of the system is repre- 
sented in Fig. |l| As already noted in in this 
case it turns out to be convenient to choose the 
following gauge condition: 

A^{xo,0) = 0, x^A,{xo,^) = 0. (4) 

Notice that this gauge choice is possible since the 
formalism is completely gauge invariant. Within 
this gauge it is possible to express the gauge field 
in terms of the field strength tensor, 

A^{x)^ I daa"-^''^ x''Fkf,{xo,ax), 
Jo 

where n(0) — and n{i) = 1. Moreover the only 
non-vanishing contribution to the Wilson loop is 
given by the quark paths connecting x with y, 
and we have 

W{T; A) ^TiP expUg J\z''Af,{z)\ . (5) 



= i-T/2, 



-T/2,6) 



Figure 1. The Wilson loop in the static limit of 
the heavy quark. 

As shown in (l^Jlltl in order to evaluate Eq. @ 
we need to know the Wilson loop average over 
the gauge fields. We evaluate it via the cumulant 
expansion described in j^. Keeping only bilocal 
cumulants we obtain: 

{W{T,A)) 

= cxp I - ^ dx'^ J' dy'^D^, (x' , y') } , 



X / d/3 (Ffc^(a;0,ax)Fj,(yO,/3y)), (6) 
Jo 

Assumption (|^) corresponds to the so-called sto- 
chastic vacuum model. Inserting expression (|^) in 
Eq. (j^) and expanding the exponential we obtain 
the following expression for the propagator So of 
the quark (which is Ginv "projected" on the first 
fermion line): 



Sd = So + Sq K Sq + Sq K So K Sq 



(7) 



So is the free fermion propagator. Taking into 
account only the first planar graph (since we 
are interested only in contributions proportional 
to the gluon condensate), we have K{y',x') = 
j'^ So{y' ,x')j'^D^^{x' ,y'). A graphical represen- 
tation of K is given in Fig. ^. Eq. (0) can be 
written in closed form as Sd = Sq + SqKSd (or 
in terms of the wave-function, {j) — m — iK)^ = 0; 
m = mi). Therefore, K can be interpreted as 
the interaction kernel of the Dirac equation as- 
sociated with the motion of a quark in the field 
generated by an infinitely heavy antiquark. 

Assuming that the correlator {F^x{x)F^p{y)) 
depends only on the difference between the coor- 
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Figure 2. The interaction kernel K. 

dinates, we define: 

(Ffe^(x",ax)F,,(y",/3y)) 

= /fcp/i.(x° -y°,Q;x-/?y). 

With this assumption K can be written in mo- 
mentum space as: 

K{q,p) 

dr / daa"('^) 

-oo JQ 

x-i" {e'(-T)A+(t)7°e-*(^'°-^')^ 

- 6'(T)A_(t)7°e-'(f''+^')^}7^ 
x/fcp/,.(T, (a - /3)r), (8) 



where t = (/3p - aq)/(/3 - a), Et = 



and A±(t) 



± ( m-t • 7)70 



Equation is our basic expression. It con- 
tains the perturbative interaction up to order 
and the non-perturbative one carried by a single 
insertion of a second order cumulant. From now 
on we want to focus our attention only on the 
purely non-perturbative interaction. The Lorentz 
structure of the non-perturbative relevant part of 
JfiXyp is 

It^uM = ]lc^^^(.9,.gAp - 9,p9Xu)D{x^) (9) 

where (F^(0)) is the giuon condensate, Ic the 
identity matrix of SU(3) and Z? is a non-per- 
turbative form factor normalized to unit at the 
origin. Lattice simulations have shown that D 
falls off exponentially (in Euclidean space-time) 
at long distances with a correlation length a"^ ^ 



(1 GeV)-i [|l3[. As shown in |] this behaviour 
of D is sufficient to give confinement at least in 
some kinematic regions. 

We notice here that the most general form of 
f'p.Xvp (^) could also contain a term of the type 

{F\0)) 1 



2ANr 



+d\{xpg^^ - x^g^p) jDi{x^), (10) 

where Di is another unknown non-perturbative 
form factor. Since this is the Lorentz structure 
of the perturbative part of /, it is not surprising 
to discover that, substituting Eq. ( |l0|) inside the 
Wilson loop in the large time separation limit, the 
terms depending on Di give rise to the kernel 

which is of the Coulomb type. The treatment of 
this kernel is trivial and gives contribution of the 
perturbative type. Therefore in this work we will 
not consider terms containing Di. 

In what follows we study expression (||) for dif- 
ferent choices of the parameters which are the cor- 
relation length a, the mass m, the binding energy 
(p" — m) and the momentum transfer (p — q) . 

A. Heavy quark potential case (m > a > p^—m) 
If we assume a to be bigger than the binding en- 
ergy (p° — m) and smaller than the mass m of 
the quark, since a ~ 1 GeV, the quark turns 
out to be sufficiently heavy to be considered non- 
relativistic. In order to obtain the 1/m? potential 
we can neglect the "negative energy states" con- 
tributions to (H) by writing 

X(<Z,p)^-g2(2^)%"-9°) 

/•oo pi d d 
X dr daa''(^^ / d/3 — — 

Jo Jo Jo op oq 



q)r 



X7''A+(t)7V/fe^.(T,(«-/3)r). (11) 

Now, inserting Eq. (^) and by means of the usual 
reduction techniques, we obtain up to order 1/rn? 
the static and spin dependent potential 

,(F2(0)) 



V{r) = g' 



2ANc 



dr / d\ 
'0 
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x(r - X)D{t^ - A^) 

dr / dX 

oc Jo 



a -LI ,{F^m 
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x[2--l]Dir^-X'). (12) 



4m^ r' 



This result agrees with the one body hmit of 
the potential given in [^,1). In particular for 
r 00 identifying the string tension a = 



mesons with a strange quark (like Dg and Bg). 
In this case the light quark mass is smaller than 
a: vis 200 MeV < 1 GeV. Actually the case 
a > m has to be considered as the only realistic 
one concerning heavy-light mesons. Under this 
condition either the exponent (p*^ — Et) as well 
as (p° + Et) can be neglected with respect to a. 
Therefore we have: 
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^ drj^ dXD{r^ - A^) we obtain ^ f' ^^^ni,) f 



the well-known Eichten and Feinberg result ||1|, 
a - Li a 



+00 

X 

'0 JQ Jo 



Op'^ Oq'- 



V{r) =ar-C 



(13) 



where C = g' 



24iV, 



dT I dXXD{T^ 
-00 Jo 



A ). We observe that the Lorentz structure which 
gives origin to the negative sign in front of the 
spin-orbit potential in ( p^ is in our case not sim- 
ply a scalar {K ar). 

B. Sum rules case {a < p^ — m, a < m) Let us 
consider now the case in which the binding energy 
of the quark is bigger than the correlation length, 
which can be considered zero respect to all the 
scales of the problem. We have 

K{q,p)^^g\2n)6{p'-q')l^I^ 



JO JO dpf" dq^ 

x{YSo{p)l^{2T:f5^{V-^))- (14) 

In particular from Eq. ( p^ we obtain the well- 
known leading contribution to the heavy quark 
condensate ||l4|| : 



{QQ) 



d'^p f d'^q 



(27r)4 J (27r)4 
1 {aF'm 
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Ti{So{q)K{q,p)So{p)} 



(15) 



C. Light quark case (a > m) Since we have 
reproduced the known results concerning heavy 
quarks, Eq. (||) should maintain some physi- 
cal meaning also when considering heavy-light 



X / dVe'(P-'i)-V(A+(t)- A_(t))7°7^ 



xtt(r,(a-/3)r). (16) 

We observe that in the zero mass limit this ex- 
pression gives a chirally symmetric interaction 
(while a purely scalar interaction breaks chiral 
symmetry at any mass scale) . This means on one 
side that our interaction keeps the main feature of 
QCD i. e. in the zero mass limit chiral symmetry 
is broken only spontaneously. On the other side 
this seems to suggest that for very light quarks 
the projectors A+ and A_ which appear in ( p^ ) 
should be taken from the chiral broken solution 
of the corresponding Dyson-Schwinger equation. 

3. CONCLUSIONS 

In the literature, also recently, a Dirac equation 
with scalar confining kernel (i. e. ~ or) has 
been used in order to evaluate non-recoil contri- 
butions to the heavy- light meson spectrum |^-|^ . 
The main argument in favor of this type of ker- 
nel is the nature of the spin-orbit potential for 
heavy quarks. This turns out to have a long- 
range vanishing magnetic contribution (according 
to the Buchmiiller picture of confinement) and is 
completely described by the Thomas precession 
term. This situation is compatible with a scalar 
confining kernel. However, assuming more sophis- 
ticated confinement models with a bigger sen- 
sitivity to the intermediate distance region, the 
spin-orbit interaction has no more such a simple 
behaviour. In particular non zero corrections to 
the magnetic spin-orbit potential show up. More- 
over, the velocity dependent sector of the poten- 
tial seems not to be compatible with a scalar 
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kernel (we refer the reader to g for an exhaus- 
tive discussion). Therefore also from the point 
of view of the potential theory there are strong 
indications that the Lorentz structure of the con- 
fining kernel should be more complicate that a 
simple scalar one. This emerges also in our ap- 
proach. The kernel (H) follows simply from the 
assumption on the gauge fields dynamics given by 
Eq. and by taking only one non-perturbative 
gluon insertion on the quark fermion line. When 
performing the potential reduction of this ker- 
nel in the heavy quark case {A) we obtain ex- 
actly the expected static and spin-dependent po- 
tentials. Therefore our conclusion is that there 
exists at least one non scalar kernel which repro- 
duces for heavy quark not only the Eichten and 
Feinberg potentials in the long distances limit, 
but also the entire stochastic vacuum model spin- 
dependent potential. Moreover when considering 
a, the inverse of the correlation length, small with 
respect to all the energy scales (case B), the ker- 
nel (H) gives back the leading heavy quark sum 
rules results. It is possible to try to extend the 
range of applicability of Eq. (^ to more realis- 
tic cases, like Dg and Bg mesons where the light 
quark mass is smaller than the characteristic cor- 
relation length of the two point cumulant (case 
C). The relevant part of the kernel is also in this 
case not a simply scalar one. We mention that 
a similar picture of the Lorentz structure of the 
confining Dirac kernel emerges also in the differ- 
ent approach of . 

An attempt to extend the present approach to 
the two-body case is given in . The equivalent 
graphs of Fig. |2| seem to play a crucial role (in 
the two-body case such kind of graph exists for 
any fermion line and as exchange graph). Never- 
theless these graphs are not sufficient in order to 
provide a complete relativistic description of the 
two-body system. The main difficulty is that in 
this case it does not exist a gauge like (||) which 
automatically cancels the contributions coming 
from the end-point strings. These contributions 
are necessary in order to restore gauge invariance. 
From this point of view the situation seems to 
be quite more complicate than in the heavy-light 
system which remains the most natural context 
to test the formalism. 
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